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Abstract 

The subject report discusses a general three-dimensional ciystal-lattice 
diffraction analysis for uniform planar antenna arrays and treats, in particular, 
arrays whose lattice phase gradients are determined by an image plane and 
simultaneous beam-forming matrices composed only of coupler hybrid junctions 
and fixed phase shifters. A large selection of array problems is discussed in a 
unified manner by means of lattice notation and the scattering matrix concept. 
Spatial coordinates of antenna beam maxima and nulls are located directly by 
means of expressions resembling M. von Laue's crystal diffraction equations. 
The subject of array gain is related to the surface integral / E 2 (<9,<jp) sin <9d#d<p. 
Topics such as grating-lobe formation, cosine-order beams, and beams in the 
"invisible" region are treated with relative ease. Weighted addition of beams 
is also considered with a view toward side-lobe level reduction for simultaneous 
beam-forming matrices, continuous beam steering, and dual-plane phased array 
monopulse angle tracking. 

This report represents an attempt to formulate the uniform planar phased 
array problem and the multiple-beam matrix concept for purposes of machine 
calculation. The scattering matrix notation introduced herein provides an orga- 
nized analysis of beam-forming circuits. It is possible to convert the scattering 
matrix to a transfer matrix which admits chain multiplication for tandem cir- 
cuits, thereby providing a method for handling circuits of greater complexity. 
The report is by no means complete and treats only a few simple scattering 
circuits. It is, however, a point of departure for a more general investigation 
of scattering junctions, circuit topology, and inertialess or electronic antenna 
beam scanning by a purely passive technique. 
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AN ANALYSIS OF 

SIMULTANEOUS BEAM FORMATION BY MEANS OF PASSIVE 
SCATTERING MATRICES AND UNIFORM PLANAR ARRAYS 

by 

R. F. Schmidt 


INTRODUCTION 

The subject of simultaneous beam formation from a single aperture has 
received considerable attention due to (1) the inherent advantages accruing from 
redundant use of the same antenna structure, (2) the inertialess scanning feature 
which eliminates physical rotation of large-aperture structures, and (3) the 
track-while-scan capability for a multiple target situation. Although high data 
rate systems are frequently associated with ballistic-missile saturation attacks 
and related tactical problems, it is evident that simultaneous beam formation 
offers numerous advantages and additional degrees of freedom for satellite 
tracking and telemetry data reception. Some of the constraints imposed in the 
tactical radar case can be removed for satellite and spacecraft problems. For 
example, the target analogue ordinarily carries a beacon or transponder and 
the array complexity introduced by short-pulse utilization is eliminated. The 
approach to multiple satellite tracking outlined is considerably simplified. 
Transient array analysis is considered to be beyond the scope of this report. 

The principal objective of this report is to present the general three- 
dimensional crystal diffraction notation and apply it to studies of uniform planar 
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antenna arrays placed over perfectly conducting ground planes, a particular set 
of three-dimensional problems. Actually, numerous uniform array problems 
can be solved in an organized manner by means of the general diffraction analysis, 
and several advantages accrue from this rather formal approach. Some of these 
are the computing and plotting of far-field radiation patterns by machine methods, 
and the determination of total radiation by means of surface integrals to deter- 
mine precise directive gain values for an array. The general three-dimensional 
analysis presented in this report draws heavily upon sources in the classical 
literature. Various assumptions and idealizations under which the analysis 
presented is valid are discussed in detail. 

The analysis is shown to provide a high order of flexibility for studying 
arbitrarily large arrays, and admits variation of such parameters as the array 
phase gradients, the element separation, the free-space element pattern, and 
the height of the source element above the reflector. Calculation of the angular 
coordinates of the principal beam maxima is achieved without mapping the com- 
plete radiation patterns, although the detailed three-dimensional Fraunhofer 
patterns can be obtained directly from simple source and lattice factors. The 
beam principal maxima and minima are located by means of von Laue type 
diffraction equations which, when satisfied, lead to simultaneous maximization 
of certain lattice factors. Using the same general three-dimensional crystal 
notation, the origin and significance of grating lobes is discussed in relation 
to so-called spurious (non-visible) beam conditions which can arise for planar 
arrays. The subject of array phase centers is introduced and an expression is 
derived for the requisite phase corrections when two or more beams are combined 
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by redundant utilization of a linear or planar aperture to form simple sum, 
difference, and cosine-order beams. 

A certain amount of material pertaining to hybrid junctions and scattering 
matrices is introduced subsequent to presenting the crystal lattice notation to 
make the application of the lattice notation specific. The array phase gradients 
are determined by particular forms of scattering matrices which, to a certain 
extent, predetermines the solutions for the far-field radiation patterns. The 
scattering matrix concept is shown to provide a convenient means of displaying 
linear and planar array gradients which greatly facilitates multiport circuit 
analysis. Basic hybrid building-blocks of the corporate feeds and the complex 
feeds themselves are described in terms of mathematical matrices whose 
individual elements represent reflection coefficients and transmission coefficients 
of definite amplitude and phase. A linear (16 x 1) scattering matrix is presented 
in detail. The corresponding simple sum beam (first-order cosine beam), the 
difference beam, and the aperture illumination (current distribution) are then 
formed to illustrate the application of the notation. 

DIFFRACTION MODEL 

Since there exists a remarkably close analogy between x-ray diffraction 
analysis and antenna array theory, the methods developed for crystal structure 
analysis are in large part applicable to antenna design. The x-ray diffraction 
problem, simply stated, is one of determining the location of the atomic dipoles 

m 

in a crystal -lattice structure from diffraction spectra consisting of intensity 

y information only. The antenna array problem is one of determining a radiation 

pattern corresponding to particular lattice spacings and prescribed phase 
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conditions among the source elements. References pertaining to the crystal- 
lattice approach, as adapted to antenna analysis, can be found in the textbooks 
of Stratton 1 and Silver. 2 In order to facilitate the discussion the notation of 
Stratton will be used throughout this report wherever the lattice analysis appears. 

Certain assumptions are fundamental to the determination of far-field radia- 
tion patterns by means of diffraction analysis. Among these is the assumption 
that the constitutive parameters are independent of field strength, which ensures 
the linearity of the constitutive equations and Maxwell's equations. That is, 
if the conductivity a, electric inductive capacity e , and magnetic inductive capa- 
city p. are independent of field strength, then the superposition of fields is admis- 
sible. A set of electric and magnetic field vectors E x and H p due to charge and 
current source functions p x and J , satisfies the field equations. This set can be 
summed with a similar set E 2 and H 2 , due to p 2 and J 2 , which also satisfies the 
field equations. Then E = E x + E 2 , H = H x + H 2 , p = p x + p 2 and J = J x + J 2 is a 
superposition which is consistent with Maxwell's equations and which describes 
a possible electromagnetic field. The subject of antenna gain is closely associated 
with tliis notion, as will appear later. 

Another assumption ordinarily taken into the diffraction analysis is that the 
radial distance r from a source element Sj to the far-field point P is given by 

r = R - t. • R 0 . (1) 

The vectors from coordinate origin to radiator, radiator to field point, and 
origin to field point are given by 7. , T and R respectively as shown in Figure 1, 
and R 0 = R/|r| , a unit vector. This notation will be consistent throughout the discussion. 
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The inference of Eq. (1) is that r x R = 

0, or 7 and R are parallel for the field point. 

Superposition of the component fields con- 
tributed by sources (1) with separate current 
distributions and (2) occupying a finite region 
of space leads to quasi-point-source fields. 

The far fields obtained in this manner are 
not point-source fields since the equiphase 
surfaces are not the family of spheres of constant radius R, in general, even if 
the individual sources are true point-source radiators. A rather lengthy argu- 
ment is required to develop the expressions for the far-field resulting from a 
multiplicity of sources if Maxwell’s equations are taken as the starting point. 
The text by Silver 3 treats this problem rigorously and in considerable detail. 

It is sufficient here to state that the far-zone fields at a point P have the form 


Z 



Figure 1 -Origin of coordinates, 
source, and field points. 
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under 


r 


= R -7. R 

j 


o ’ 


where the new terms are: 

T = the magnetic current density; a formalism taken to be zero, generally, 

J m 

but retained here for symmetry purposes. 
oo = 2n f 
f = frequency 



k = wavelength 
i = v~ 

The volume integrals are independent of r, ensuring that RE p and RH p remain 
finite as the radius of the sphere of observation approaches infinity (R - oo). 

Since E and H satisfy the radiation conditions 

R-"!» R [(*o xh )* (j)‘ /!g ] = 0 (4) 

and 

«-- R [0 ,/, <s*“ E >- R ] =° (5) 

the radiative components of E and H are mutually perpendicular and both normal 
to R 0 . 
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At this point of the development, it is sufficient to consider only the E fields 
since radiation has now been established, and it is helpful to resolve the E field 
of Eq. (2) into <9 and 0 components in a spherical coordinate system. Then, if 
i e , are unit vectors , 


E p = E e l e + 


(6) 


where 


1 ^ „-ikR 

4^ 


f r~ . / 

e \i/2 _ . ~ 

i k r . *R n 



'jjj Jm 

e J 0 dv 


'v 




Fj ( 0 , 0) 



(7) 


and 


E <* = E p' i * = 


477R 


I 

m 

% T / 

c\l/2 - 

e~ lkR 

V 

} ‘ l * " ( 

fl) J"i’ l 0 


dv 


„-ikR 

4ttR 


F, (&, <t>) 


( 8 ) 


since 


J m * R 0 • h = L • K 


and 


L x Rn • L s - L ■ C 


'0 V = ■ Jm e 


(9) 


( 10 ) 


The diffraction analysis for a planar array can now be accomplished by consid- 
ering the cross-polarized components Eg, E^ independently, therefore, the 
crystal-diffraction analysis will be written for linear polarization only. Com- 
posite fields can then be constructed as required by summing the resultant Eg and 
E^ fields radiated from the array with due regard for their phase relationship. 
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The amount of complexity introduced into an antenna array analysis depends 
upon the assumptions admitted for the ensemble of radiators or scatterers. In 
the case of x-ray crystallography, the scattering cross-section is a small frac- 
tion of the area irradiated by the incident wave and the usual assumption is that 
the scattering from any one center is independent of the presence of other scat- 
terers. 4 Multiple scattering between elements is neglected. This assumption is 
carried into the present analysis and each element of the array will be consid- 
ered identical and invariant with scan angle. Each element is also taken to be an 
interior element of the array, to the exclusion of all edge effects. An equivalent 
statement is that the array is assumed to be of such an extent that the ratio of 
peripheral to interior elements is negligibly small. Any interaction which leads 
to a new set of identical element currents IV will not invalidate the analysis 
since the original set of currents I is simply replaced by I V which permits 
a reapplication of the notation. Problems pertaining to power conservation, 
energy storage in the array, absolute gain, or absolute field intensity require 
additional attention. 

LATTICE NOTATION 

A planar array of sources exhibiting complex phase patterns, situated over 
a perfectly conducting ground-plane is a particularly interesting array which can 
be considered as a three-dimensional lattice of source elements and is one which 
brings out most of the salient features of the analysis. It will be assumed that 
the array is uniform in the sense that (1) all of the real and image elements are 
parallel to one another and carry currents of identical magnitude, (2) the spacing 
of the elements is uniform but not necessarily equal on each of three orthogonal 
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coordinate axes, and (3) the phase gradient on each axis is linear, but not neces- 
sarily equal for the three principal directions , a^, as shown in Figure 2. 
The result of Eq. (7) can be restated, with some slight modification, to admit 
complex array elements and to provide for relative phase delays between the 
array phase gradients. It is convenient to undertake this by noting that the j th 
oscillator in the lattice is located by a vector 

r j = i i a i + h a 2 + j 3 a 3 

as shown in Figure 2; where j , j 2 , j are integers including zero. It follows 
that the phase of the j th oscillator with respect to the phase of the oscillator at 
the origin is fi . , 


Pj = (j i a i + i 2 a 2 + i 3 a 3> + ( S i + S 2 + 5 3>> ( n ) 

allowing for phase shifts Sj , § 2 , in the network corresponding to the forma- 
tion of the gradients a l ,a 2 ,a 3 respectively. Following Stratton's 5 lattice nota- 
tion, the radiation field of the j th oscillator is 

ikR” i^t ~ ikRQ*rj ~ i/3. 

E. = - i 60 I . F n (0) 2 (12) 

J Oj 0 V > p 

assuming a half-wave oscillator oriented as in Figure 2, where R Q is a unit 
vector, R Q -R/ |r| , in the direction R and is consistent with Figure 1. The 
composite or resultant field for the entire array then becomes 


E = - i 60 F 0 (6>) 


e i kR - ico t 

R 




(13) 


If the simple half-wave oscillator is replaced by a complex source element, and 
if one of the E 0 , E^ fields is considered at a time, then Eq. (13) can be written as 
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Figure 2— Uniform lattice structure in 
three dimensions defined by orthogonal 
base vectors a^, a^. 


E = - i 60 F'(0) e _i «( 9 ^) 


— F I e " i(k6 °‘ 7 i + 

1 Oj 


where F' (0) represents the amplitude envelope of the free-space field-intensity 
pattern and g(<9, 4>) is a function which defines the phase-pattern of the com- 
plex source element. 

Let the number of radiators in any line parallel to the x, y, z axes be n x , 
r> 2 , n 3 respectively, so that the lattice is completely filled, and let the spacings 
along the coordinate directions be a t = |aj , a 2 = |a’ 2 | , a 3 = |a 3 | . Unit vectors 
are then given by 


a a, 

1 = — • 1 =— , 1 
x a 1 y a 2 2 


- - l a i \ l a 2 \ / a, 

R n 'r . - sin 6 cos 0 — + sin 0 sin 0 — + cos 6M — 

J I \ a i/ \ a 2 J \ a 3 


a i\ . /M 
J l a l [-J +} l a l {— 2 ) + 


sin 6 cos 0 + j 9 a 9 sin 6 sin <$> + j . a.. cos 0 , 


Complex lattice factors can now be defined as follows to account for the total 
field generated by the array of oscillators. 


f 


1 



— ijj (k sin#cos<£ + o.j) 


( 16 ) 


f 


2 



— i j 2(k a 2 sin 0 sin<£ + a 2) 


(17) 


f 


3 



-i j , (k 


cos $ + a 3 ) 


(18) 


Each lattice factor above has exponentials whose arguments depend upon a spa- 
tial term (a projection) and temporal term (a phase gradient). All diffraction 
phenomena discussed in this report are based on these fundamental expressions. 


Let 


)'l = k 3j sin 6 cos 4> + 
y 2 = k a 2 s i n 6 s i n + « 2 

y 3 = k a 3 cos 6 + » 3 

so that three geometric progressions can be identified as follows. 


f 



(s = 1, 2, 3). (n s > 1). 


(19) 

( 20 ) 
( 21 ) 


( 22 ) 


The series does not converge in general, but for a finite number of terms, 


-i y -i 27 

f = e° + e s + e s + 


+ e 


-i(n - 1 ) y s 
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and 


— i 7 -iy -i 2y -i(n ) y , 

(f e s ) = e s + e s + + e • " 


therefore, 


-i y — i n 7 

f - f e s = 1 - e s s 


so that a closed-form expression can be employed for the uniform array such 
that 


, “in 7 
1 - e s s 


-in 7/2 SI 
- S S 


1 - e 


-i y 


y / 2 


■f£) 


— 1 ( n -1)7 /2 

e s sin (n y/2) 


sin (y /2) 


(23) 

(s = 1, 2,3) 


Equation (14) can now be rewritten as 


E = - i 60 I 0 F;(d?) elile i(S i +5 2 + 8 3> ^ f 2 f 3 e ikR - ic<Jt (24) 

R 

This equation is sufficient to map all patterns for three-dimensional uniform 
arrays. 


The collection i 1 , f 2 , f 3 constitutes the lattice factors of the array and the 
magnitude of these factors (eq. 16, 17, 18) is F , F , F 3 respectively. Since the 
lattice factors are associated with the amplitude or field intensity patterns in 
subsequent work, it is useful to determine the maximum value of each F 2 , F 2 , F 3 
and the least upper bound (sup remum) of F X F 2 F 3 . 


F 



sin (n y 2 

v s s 

sin (y s /2) 


c -.(n s -l) 7/2) 


(s 


which is indeterminate for 


1, 2, 3) 


(25) 


7/2 = 0, ±77, ±277, 


(h s rr) 
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where h s = any integer. 

lim sin(n s y/2) lim n/2 cos(n s y s /2) 

y s /2~hn sin(y s /2) “ 7 s /2 - hw l/2cos(y s /2) ' (26) 

The maxima of the lattice factors are therefore given by rij , n 2 , n 3 correspond- 
ing to fj , f 2 , f 3 respectively, and the least upper bound for the summation of 
component fields occurs when ^ F 2 F 3 = n 3 n 2 n 3 . This condition, simultaneous 
maxima for F x , F 2 , F 3 is not a prerequisite for obtaining a principal beam 
maximum. That is, E (principal beam) < E (supremum) in general. 

The function sin(n s y s /2)/sin(y s /2) can be positive or negative, depending 
on the arguments (n s y s /2) and (y s /2) , and is of considerable interest in gen- 
eral array studies. A plot of the numerator, denominator, quotient (field pattern) 
and quotient-squared (power pattern) is shown in Figure 3. Particular significance 
attaches to the sign of this function, or more precisely, the sign of the product of 
three such functions without the exponentials e l(r,s 1)r « /2 ( s = 1,2,3). It is 
shown later in this report that the summation of simultaneously formed beams 
ultimately depends upon the sign of the triple product of lattice factors for each 
beam taken into a sum via the field superposition theorem. A detailed discus- 
sion of this, including the removal of the exponential terms e * (ns t e ~*^ 

arising from the lattice factors and the source elements, respectively, is de- 
ferred to the discussion of scattering matrix gradients and beam summation. 

ARRAY GAIN 

An approximation to the maximum array gain can be made from Eq. (24), 
based on the principle of field superposition, if it is assumed that each of the 
n 2 n 3 real (physically energized) sources dissipates as much energy when in 
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array over a perfect ground plane (n 1 2) as when placed in a free-space en- 

vironment. For definiteness half-wave dipoles are selected as source elements. 


The general expression for antenna gain 6 is 


G = 


P('\ 4 ) 



P(r, 0) sin tkl&dcfr 


(27) 


relative to a hypothetical isotropic radiator, from which the power gain of a half- 
wave dipole in free space is 


G 


X/2 “ 


1 


4 77 



sin 6 J d(^d0 


(28) 
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where 


F Q = cos cos 6 


sin 6 



(F q ) 2 sin^d^d^ = 8 = 


477 

1.65 


( 29 ) 


is the total energy radiated from a single \/2 dipole. The validity of (29) fol- 
lows from 

1 2 

G X/ „ . — 1- r = 1.65 

A./2( max . ) j \ 

477 \ 1 . 65 / 


which implies 2.17 db gain, a well-known value. The total energy S T from n 2 n 3 
dipoles of equal illumination is 




(30) 


which represents the total energy from the array by assumption, 
gain is 


477 


^ irk) 


1.65 (F„F,F 2 F 3 ) ! 


11 2 n 3 


Then the array 


The least upper bound (supremum) for all possible values of G ft is obtained for 
F q = 1, Fj =2, F 2 = n 2 , F 3 = n 3 - 

C „,„p> = < 4 ' < 165 ' "2"3 


G., ... = 6.02 db. + 2.18 db. 

A ( s up )d b 


+ 10 log 10 n 2 n 3 db. 


By way of example, for a square array of 256 elements, n 2 n 3 - 256 and with 
a 2 = a 3 - k/2 say, 

G (sup) db = 8. 19 + 10 log in 256 = 32db. 

16x16 10 
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Actually, the total power radiated from a collection of oscillators in array 
carrying an illumination I Qj is not, in general, equal to the totality of real 
oscillators times the power radiated from one such oscillator with the same 
illumination when the latter is considered separately. 7 Since Maxwell's equa- 
tions are linear the superposition of fields is always admissible, but (1) addition 
of the Poynting vector (P = EXH)or (2) addition of source element power, pertain 
to nonlinear combinations and are incorrect solutions for total power except in 
a few very special cases. It has been pointed out by Stone 7 and others that a 
pair of atoms radiate independently for a phase gradient a 3 - v/2 , separation 
a 3 arbitrary and behave almost independently for a = 0, a = tt if a 3 > 0.7\. 

These special cases have microwave analogues in a pair of identical linear and 
collinear antennas. 

In general, the comparison of diffraction patterns should be predicated upon 
equal input power. This can be accomplished by (1) deriving the shape E-(F Q F 1 F 2 F 3 ) 
of the patterns by means of the diffraction analysis, (2) forming the Poynting 
vector (P - EXH) , (3) determining the total radiation (P .= j Ps inPdfc-'dc/; = K J E 2 ds) 

S const. S 

and (4) normalizing the patterns with respect to total power in each (■' r K J 

s 

(/C^ E 2 ) 2 dS , where Pj - C Q P relates the total power in each pattern). The field 
patterns can then be weighted by r C Q which allows pattern comparison on an ab- 
solute basis. The subject of directive gain is presently being investigated with 
the assistance of an IBM 7090 computer to approximate the integral of Eq. (27) 
for various arrays under different conditions. The effect of element spacing, 
aperture gradients, and beam addition on antenna gain will be determined on an 
individual basis in subsequent reports for several array configurations. 
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VON LAUE EQUATIONS 

An appreciable amount of information can now be derived from the crystal 
lattice notation including (1) position of principal beam maxima, (2) grating lobe 
information, and (3) spurious beam information. Returning to Eq. (26), the in- 
determinacy leading to the realization of the maximum value for each lattice 
factor is identified by 


s in y s /2 = 0 (31) 

which occurs when y s /2 = h s 7 r if h s equals an integer (plus, minus, or zero). 
Then, from Eq. (19), (20), (21); 


y 1 = kaj sin 6 cos <f> + = 27rh 1 


or 


y 2 = ka 2 sin sin 0 + a 2 - 27rh 2 
y 3 = ka 3 cos 0 + a 3 = 27 /h 3 


s in 6 cos 4- 




(32) 


sin 0 sin ^ t (h 2 1 J A u 2 (33) 

\ 277 / a 2 


cos 0 = fh, — ] — - n, . (34) 

\ 3 2tt) a 3 

Equations (32), (33), (34) are the array analogues to the equations obtained 
by Max von Laue 8 for crystal diffraction maxima. These conditions are not 
necessarily satisfied simultaneously by any choice of 0 and </> coordinate angles 
for arbitrarily chosen hj, h 2 , h 3 . The choice of the lattice spacing and the array 
gradients determines the possibility of simultaneous solution of these equations. 
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A few simple examples illustrate the usefulness of the von Laue type equa- 
tions. An array of elements which forms images by means of a perfectly con- 
ducting ground plane introduces a specific n s and a specific a s . If the plane of 


the array is taken to be the y-z plane, then n 3 = 2 and - tt for source ele- 
ments oriented parallel to the y-z plane due to boundary conditions at the ground 
plane. Simultaneous solution of Eq. (32), Eq. (33), and Eq. (34) yields 


2 2 2 
u i + u 2 + U 3 = 


h, - iU 

W a. 




h 2 - -i)i. 

2nl a. 


-i2 r 


h 3 - — ) — 

2ttJ a 3 _ 


(35) 


If aj = a 2 = a 3 - k/2 is arbitrarily selected for the array with source ele- 
ments over the ground plane, then 0 < 0 < rr , - rr/2 5 0 5 + 7 t/2 for physically 
realizable beams and 

1 = (2h, - 1)’ t t ^2h, - ^ . (36) 

The first term of Eq. (36) is minimized for h x - + 1 or h x = 0, the latter value 
being non-realizable for the array and ground-plane combination here since the 
beam forms in the physically inaccessible half-space. Then 

1 ^ 1 + ^2h 

The two parenthetical terms must vanish if Eq. (37) is to be satisfied. If a 2 = 
a 3 = 0 is imposed on the array there exists h 2 = h 3 = 0 such that equality holds. 

If a 2 = a 3 = 7 t 16 is imposed on the array there does not exist any such combina- 
tion of h 2 and h 3 and inequality holds for Eq. (37) 




(37) 


Since « 2 and a 3 are the gradients ordinarily associated with beam scanning, 
for the assumed array orientation, it is desirable that these parameters take on 
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many values besides 0 , ±77, ±27 t , etc. It is apparent that Eq. (35) cannot be satis- 
fied for all choices of a r , a 2 , a 3 and 0 ^, a 2 , a 3 . The factor F Q f in Eq. (24) can 
be regarded as the pattern resulting from the free-space source element F' and 
the influence of the ground plane fj , so that the array problem reduces to a study 
of finding simultaneous solutions to Eq. (33) and Eq. (34) only. That is, a field 
point P(6, 0) is to be located such that the maximum values F 2 = |f 2 | = n 2 

m m 

and F 3 = I f 3 j = n 3 are to be brought into coincidence at some coordinate 

m m 

point P(0, 0) and a value F x - |fj < F 1 - n 1 , will be accepted. F t < F x im- 

m m 

plies that the beam principal maximum, if it exists, will be less than the least 
upper bound (supremum) of all the possible values of E attainable under free 
choice of , a 2 , a 3 and c^, « 2 , a 3> 


Simultaneous solution of Eq. (33) and Eq. (34) yields 



s in 0 = 

(1-uj) 1 ' 3 



under - 0 < 6 < 77, 0 £ 0 5 2-n 

which provides the angular coordinates ./) if they exist. 


(38) 


(39) 


Three simple examples are given to illustrate the significance of Eq. (38) 
and Eq. (39) in relation to normal (ordinary) beams, grating lobes, and spurious 
beams. (See Figure 2.) 
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Example (I) 

Let cij =77, aj = a 2 = a 3 = X/2 and 0 < 0 < 77/2, -rr/2 < 4> < +77/2 

a 2 = a 3 = 77/1 6 is imposed. 

h 2 = h 3 = 0 are unique possibilities. 

One principal maximum forms, subject to F 0 Fj. This is an example of a 
"normal" beam. The array is over a perfectly conducting ground plane, the 
separation between all real and image oscillators is \/2, and equal gradients 
are imposed in the (y) and (z) directions. The beam is "normal" in the sense 
that it is one of set of beams ordinarily formed by means of scattering matrices 
which generate gradients a 2 = u77/n 2 and a 3 •- v77/n 3 , where u and v are odd 
integers. "Normal" or "orthogonal" beams as they are sometimes called have 
the characteristic that the principal maximum of one beam coincides in space 
with the minimum of an adjacent beam. 

Example (II) 

Let ctj = 77, a t = a 2 = a 3 = K and 0 < 6 < tt/2 - n / 2 £ <£ £ + rr/2 

a 2 - 77 / 16, a 3 - 13 7 r 16 is imposed. 

h 3 * 0, h 3 - + 1, h 2 = 0 are unique possibilities. 

Two principal maxima form, subject to F Q F r In this case a so-called 
"grating lobe" forms, but the principal maximum does not coincide with the 
grating lobe maximum in space. The array is over a perfectly conducting ground 
plane, the separation between all real and image oscillators is X and unequal 
phase gradients are imposed in the (y) and (z) directions. In this instance the 
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increased element spacing allows the formation of grating lobes, if a 2 or a 3 
are made sufficiently large. Grating lobes exhibit intensities which are equal to 
or less than the principal maxima for an array unless the source pattern F Q and 
the lattice factor F associated with it diminish the product F Q Fj F 2 F 3 as is the 
case for large steering angles. Figure 4 illustrates typical linear grating lobe 
formation after F Q F 1 is allowed to operate on the factor F 3 . F 2 = 1 for a linear 
array here. 




(b) a 3 - 15 77/16 RADIANS 


Figure 4-Principal maximum (T), and grating lobe (2), of linear array over ground-plane; 
a j = A/2, a 3 = 6 A/10, and A/ 2 dipole source elements. 
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Example (ID) 

Let aj = 77. = a 2 = a 3 = A./ 2 

a 2 ~ a 3 ~ IStt/16 0 £ 6 < 77 , -77/2 < £ +77/2 is imposed. 

A pair h 2 , h 3 does not exist. "Spurious beams" result instead of principal 
maxima, subject to F 0 F X . The array is the same as in the first example, except 
that the phase gradients are larger. In this instance there does not exist a point 
in space such that the maximum value of the product F 2 F 3 is realized. That is, 
the last two von Laue equations, Eq. (33) and Eq. (34), cannot be solved simul- 
taneously for any given direction in space. This is sometimes referred to as 
"spurious" beam formation. 

The von Laue type equations lead to a precise determination of the angular 
coordinates ( 6 , 4 >) of the beam principal maxima where F 0 F X is isotropic only, 
a condition which is approached quite closely for many choices of source ele- 
ments if the beam is not steered too far from the normal to the plane of the 
array. Ordinarily the exact field patterns are determined by evaluating 
F 0 FjF 2 F 3 , but it is useful to retain the hypothetical isotropic sources here and 
elsewhere to obtain an insight of the function of the lattice factors F 2 F 3 without 
the effects of the FqFj term. 

The lattice notation requires knowledge of the free-space pattern F 0 only, 
and the effect of the ground plane is introduced by including the term F x . As a 
specific example, the field of a single thin linear center-fed dipole in free space 
is given by 9 

160 |I | 

E - 

t> D 
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where L is the total dipole length. For the special case where L = k/2 Eq. (40) 
reduces to 



(41) 


Equation (41) corresponds to F 0 in the lattice notation and is therefore sufficient 
to describe this particular source element. It can easily be verified that 


cos f— cos 8 

\ 2 j 


s in 6 


- F F 
cri 


(42) 


leads to the same three-dimensional radiation pattern as the more complex 
expression 10 


F of,= E » = i(f) 


1/2 


I e -j [kR + (<£/2)] 


77R 


:os ^ cos 


sin 6 


cos 


(^sine-sin i-±) 


(43) 


where the following correspondences can be made between the notation of Silver 
and Stratton. 


Silver 

Stratton 

Value 

j 

i 

V' — 1 

e -jkR 

e+jkR 

- 

7 

a i 

77 

a 

a i 

- 


Equation (43) holds for the half-wave dipole a distance a 1 /2 above a perfectly 
conducting ground plane and parallel to it. The relative ease with which source 
patterns can be developed, and the fact that the size or extent of the array does 
not increase the computational effort required to evaluate the lattice factors 
makes the lattice notation especially attractive. 

*The free-space source factor F 0 for a pair of crossed dipoles in arbitrary phase relationship is 
given in the Appendix A. 
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BEAM MAXIMA, MINIMA, AND CROSSOVER LEVELS 
UNDER ORTHOGONALITY CONDITIONS 

An examination of the function sin (n s y s /2)/sin (y s /2) , as it appears in 
Eq. (23) and Eq. (24), leads to a set of equations for determining the beam nulls 
of a planar array. The procedure is similar to that for determining beam max- 
ima by von Laue type equations. It is sufficient, however, for any one of F Q , Fj , 
F 2 , F 3 , to vanish to establish a radiation pattern null or zero value for the field 
intensity. 


f = F = • 

*i l 


sin — (k Sj sin 6 cos + a^) 
sin-^- (k aj sin 0 cos ^ + a j) 


= 0 


(44) 


if 


n i 

— (k sin 6 cos ^ + o^) = 


TT 


(45) 


where ■i 1 belongs to the set of all integers (f^ez) . Then 

sin<9cosc?!>=— — ) (46) 

a i \ n i 2tt J 

with the restriction jL njhj since ^ .. nj hj leads to a maximum value for 
the lattice factor Fj . That is, Eq. (46) becomes Eq. (32). 


Similarly, 


if 


f = F 


n 2 

sin— (k a 2 s i n 0 s i n 4> + a 2 ) 


sin — (k a 2 s in 0 s in + a 2 ) 


n 2 

— (k a 2 sin 0 s in </> + a 2 ) = l 2 tt; (^ez). 


(47) 
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Then 


s in <fi - 


k 


a 2 a 2 


s in 0 l ru 2tt J 


£ n 2 ^2‘ 


(48) 


Similarly, 


if 


n„ 


sin — (k a 3 cos 1 9 + a 3 ) 

f 3 = f 3 = \ = 0 

sin — (k a 3 cos 6 + a 3 ) 


(49) 


Then 


— (k a 3 cos Q + a 3 ) = I^tt- (^ez). 


a ^ (^3 a 3\ „ , 

cos 0 = — — — ; T, p n, h 

a 3 l n 3 277/ 3 


3 “3 


(50) 


The usefulness of Eq. (46), Eq. (48), Eq. (50) can be demonstrated and some 
of the abstract character of these equations removed by an example. Assume an 
array over a ground plane such that arbitrarily chosen values are 


n i = 2, n 2 = 1 6, n 3 - 1 6 

a i - a 2 = a 3 = X/2 

U77 VTT . N 

a l= n ’ a 2=~- a 3 =— : ( U ' V ) £ ( Z odd)- 

2 3 

so that a 2 = 1-7 16, 3 tt/16, ••• , 15t 7 16 and a 3 - 177/16,3 tt 16 ••• , 15tt/ 16. 
Phase gradients given by the odd-integer rule above are physically realizable 
and lead to the concept of simultaneous orthogonal beams. A physical means 
for obtaining such gradients is presented later in this report under composite 
scattering matrices. 


From Eq. (46); -t = njhj = 0 implies Fj =F 1 (max) and hj = 0, but can be 
ignored for the half-space given by 0 - 6 - 77 and -tt/2 < £ + n/2 . = +1 
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implies that there exists an infinity of zeros corresponding to </> = +tt/2, 0 - 6 - tt 
and 4> - -77/2, 0 < 6 < n. ■i l - +2 = r^hj implies F = Fj (max) and h x = +1 which 
is in the physically accessible half-space. > 2 and ■£ <0 will not satisfy 
Eq. (46) and can be ignored. This exhausts all possible choices of and for- 
mally establishes the fact that the lattice factor Fj has only one physically re- 
alizable maximum for the assumed array. The infinity of zeros in the y-z plane 
ensures that the field intensity of the far-field pattern vanishes in the extended 
plane of the metallic boundary. The latter statement is independent of spacing 
aj in Eq. (46) and is a consequence of n = 2, a = n. 

A plot of the principle beam maxima positions for the assumed planar array 
appears as shown in Figure 5, and results from the use of the generally valid 
Eq. (38) and Eq. (39). The factor F 0 Fj can be ignored here if the source pattern 
is slowly varying with (9, 4 >) . 
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It appears from the plot of Figure 5 that the beams are equispaced in (u 2> u 3 ) 
space. This is shown by simple proof. Since the assumed array exhibits no 
grating lobes it follows that the unique values h 2 = h 3 = 0 are associated with 
simultaneous maximum values of F 2 F 3 if a beam principal maximum exists. 

The total number of principal maxima in the physically accessible half-space 
is seen to be 208 rather than n 2 n 3 = 256. Equations (33) and (34) cannot be 
solved for 48 combinations of (a 2 ,a 3 ) which lead to spurious beam formation. 
Using 

u 2 = sin £? sin = ^h 2 - , (33) 


u 3 = cos 6 





(34), (38) 


U 2 “ “ a 2 /rT ’ (51) 

u 3 = - a 3 /7r. (52) 

Since o- 2 - uir/n 2 and a 3 = v77/n 3 and n 2 = n 3 here (square array), it follows that 
"adjacent" beam principal maxima are equispaced in (u 2 ,u 3 ) space and form a 
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perfectly square grid as shown in Figure 5. "Adjacent" implies a 2 is a constant, 
and a 3 is a variable or vice versa. The coordinate system of Figure 5 is de- 
veloped graphically in Figure 6, where concentric circles in the y-z plane are 
used to map space points associated with (0 1 ,4> 1 ), ,4> A ), (< 9 x , <£ 10 ) and (& 2 ,4> x ), (0 2 , 

(d 2 >4> l0 ) or the u 2 , u 3 plane. 



It is now relatively easy to show that the nulls of adjacent beams occur at 
common angles for orthogonal beams. From Eq. (48) and a 2 = un/n 2 , 


s i n 'p = 


X 


a. sin 0 \ n 


TT-^r- ; ^2^ n 2 h 2- 


2n 


(53) 


yields a zero for the field pattern. For any two adjacent beams generated by u 
and u + 2 = u' (odd integers) the following is possible: 




U_ _ _ 2 _ 

h 2n 2 n 2 


u 

2n~ 
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Then the primed and unprimed beams have a common zero when 1' 2 = l + 1. A 
similar argument can be presented using Eq. (50). 


\ AC, a 

cos 0 = ± _L._L 


a, \ n. 


277/ ’ 


^3 ^ n 3 h 3 " 


Here a 3 — vtt/ n 3 , and v and v' = v + 2 are used. 


In the above (u ,v) e (z odd) as before. 


(54) 


It also follows that for adjacent orthogonal beams the principal maximum of 
any beam coincides with the nulls of all other beams. The proof is presented for 
Eq. (48) only since the same argument can be applied using Eq. (54). A beam 
maximum is implied by l 2 = 0 for the unprimed beam which is generated by an 
aperture phase gradient (u) rr/n 2 . A minimum is implied by i‘ 2 =1 for the ad- 
jacent or primed beam generated by a (u + 2) n/n 2 gradient. Then 


. , f\. 

s in <p = 

n 2 a 2 sin 0 

and the proof is complete. 



n 2 a 2 sin 0 


1 - 


u + 2 


The angular position <}S) at which the principal maxima of adjacent orthog- 
onal beams intersect can be found for the planar array using Eq. (24). Elimina- 
tion of factors not relevant to the crossover level determination leaves 

F o f , F 2 F 3 = k f ; f ; F 3 <55) 

where F 0 = , Fj =F' at the point of crossover only, and F 3 = F 3 if the two 

beams are adjacent in the sense that a 2 - un/r> 2 and a' = (u + 2)7r/n 2 and 
a 3 = a ' 3 . The problem then reduces to 


4 
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.n jj-^ (k a 2 sin 9 sin + a 2 )J I sin (ka 2 sin 6 sin 0 + a')J 


;in — (k a sin d sin 4> + a 2 ) 


sin— (k a 2 sin 6 sin + a 2 ) 


via Eq. (25). 


The condition for a solution to Eq. (56) is that 


(k a sin 6 sin 4> + a ) = - (k a sin 6 sin 4> + o.') 


which can be shown by the following. For small values of the argument (y /2) 
such as are encountered for beam intersections lying between the principal max- 
ima and the first null, 

sinny/2 sin n y 72 sin n, (k a, sin 0 sin ^ + a ) 

s 5 s s - f _ /[-q\ 

siny s /2 y s /2 (k a 2 sin i9 sin 0 + a 2 ) 


Let k a, sin(?sin<p=y and a 2 > a ' . Then 


f(y, a ) 


sin n 2 (x + o - 2 ) 
(X + a ) 


g(X, a') 


sin n 2 (X + «•') 

(X + <) 


which are both even functions. It follows that 


g(X, a') f(,V, a ) 


at the point of intersection of adjacent beams. But Eq. (59) and Eq. (60) are even 
functions. Therefore, Eq. (61) is satisfied if 

(X + s) r (x + ai) (62) 


(X + «- ) - _ (x + a'). 
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Eq. (62) forces the conclusion that a 2 = a' 2 which represents two beams every- 
where coincident in space, and violates the assumption that a 2 = un/ n 2 and 
a' 2 = (u + 2) 7r/n 2 for adjacent beams. From Eq. (63) it follows that 2 k a 2 sin 8 
sin 4> c - - C «2 + which is the condition stated in Eq. (57). Rewriting, 

— (cc 2 I (1.) “ (u -f 1) 

sin 4> c = = (64) 

2ka 2 sin 6 2 n 2 a 2 sin 6 

where 4> c is the value of 4> at crossover and u is the least value for the two 
beams selected. Similarly if a 2 - a' 2 , but a 3 -t a ' 3 , 


_ X. (v + 1) 

cos 6 

2 n 3 a 3 

where d is the value of 6 at crossover and v is the least value for the two 

c 

beams selected. 


(65) 


The utility of Eq. (64) and Eq. (65) is illustrated by two examples. 

Let u = -1 and u' = -1 +2; a 3 = >x 3 . Then the lesser of -1, +1 is -1 in Eq. (64). 


sin d> 

r C 


-k (u + 1) 
2n 2 a 2 sin 0 


- 0 and 0-0. 

r c 


Let v = -1 and v' = -1 +2; a 2 = a'. Then the lesser of -1, +1 is -1 in Eq. (65). 


- (V f 1 / 

cos 0 - 0 and 0 ~ v 2 . 

2n 3 a 3 

If <p c is substituted into the expression for the lattice factor F 2c ,Eq. (25) for s = 2, 
or if 6 c is used to evaluate the lattice factor F 3c , Eq. (25) for s = 3, the follow- 
ing hold. Here n 2 > > 1 or n 3 > >1. 


- 2 n 2 / 77 

(66) 

- 2 n 3 A r 

(67) 
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Assuming F 0 = F' and Fj = Fj for the adjacent beams it 

prin.max. prin. max. prin. max. prin.max. 

makes sense to state that the crossover level, referred to the main beam princi- 
pal maxima, is given by 


2c 


2m 


77m 


3c 


2n, 


7m, 


_ 2 _ 

77 


( 68 ) 


prin. max. 


prin. max. 


even if 


< 1 and 


< 2. e„ 


C 2c’ ^3 


■3c’ 




1 implicitly. 


prin. max. prin. max. 

That is, the crossover level is approximately 3.9 decibels below the beam prin- 
cipal maximum level under the special restrictions given above. 


The preceding discussion, pertaining to adjacent beam crossover levels for 
planar arrays, holds equally well for linear arrays and either Eq. (64) or Eq. (65) 
can be made applicable by proper orientation of the linear array in the coordinate 
frame. A set of E-plane radiation patterns derived for a 16 X 1 (linear) array of 
half-wave dipoles is shown in Figure 7. The difficulty of defining the beam cross- 
over levels for all steering angles in any simple manner is illustrated by the pat- 
terns with the largest beam-steering gradients. 

SIMPLE SCATTERING MATRICES 

The subject of corporate feeds is extensive and will be limited here to a brief 
discussion of passive circuits composed of hybrid coupler elements. Even this 
restricted topic can lead to numerous involved discussions, therefore, the present 
development will be carried just far enough to allow several examples of the ap- 
plication of the crystal lattice notation given earlier. Since the basic components 
of beam-forming matrices are 3-db. directional couplers and hybrid rings, the 
mathematical representation of these components is given without derivation. 
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followed by one composite scattering matrix associated with antenna beam scan- 
ning circuitry. 

Coupler and ring type "magic tees" are shown in Figure 8 together with their 
respective (coaxial) scattering matrices. 12 The convention adopted here is that 
only the angle \p , say, will be entered as a matrix element of large matrices for 
a real time delay e where the context makes the meaning clear. A scattering 
matrix element is taken as the reflection or transmission coefficient between an 
input terminal (j) and an output terminal (i) , where i j is the row-column 
designation. Both Q S.. and D S.. are symmetric matrices since S.. - S . . and 
are in fact both unitary 13 matrices since the complex conjugate (S. .)* trans- 
posed, or (S*.)*, times the original scattering matrix gives the identity 
matrix I , where I is a diagonal matrix with unit elements. 


z° 

% 

Zo 


Zo 


X 

V2 

X 

u 

Zo Zo 

/4 


Zo 



72 



X /4 



0 0 | 1 j 


0 0 1 1 -1 

0 0 J J 1 

_ j 

0 0 1 l 1 

1 ilo 0 

° s i \ 

\ 

, 

1 1 j 0 0 

i 1 ! 0 0 


-1 1 t 0 0 


COAXIAL COUPLER CIRCUIT COAXIAL RING CIRCUIT 

Figure 8-Magic T coaxial configurations and the 
associated scattering matrices. 
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Only the coupler type circuit will be employed in the remainder of this re- 
port, however, the equivalence of the ring and coupler types is formally estab- 
lished by the following transformation. The terminal (reference) planes of the 
four lines attached to the coupler circuit can be displaced by the transformation 


□ 



PS P 

□ 1 j 


(69) 


where P is a diagonal matrix with elements P k k = e‘ )l ^ k corresponding to the 
k th line into which a phase-shifting element is introduced. The choice of the 
i p k > 0, above is consistent with the use of a negative angle for a real time delay 
in the writing of the S.. previously. If the two scattering matrices are equivalent, 


o 


S. . 

* j 


□ 


S'. . 

i 1 


PSP 

□ 1 J 


(70) 


After premultiplication and postmultiplication of Q S.. by P the Q S'. is com- 
pared to 0 S. on an element -by-element basis, yielding the following equalities. 


j t e - j(0 i + ^'3 ) _ e -i<* 2 +0 4 ) ( 71 ) 

_ i ! e - j(v > + ^ ) (72) 

- 1 — e ~ j (v -’ + ' 4) (73) 

Therefore = 0, 0 2 = 7T /2, <P 3 = 3-/2, = rr , 0. > o are solutions. Place- 

ment of corresponding delay lines in lines (1), (2), (3), (4) of the coupler cir- 
cuit is required to make „S and S indistinguishable from one another. 

O ,j □ ij 


COMPOSITE SCATTERING MATRICES 

The synthesis of composite scattering matrices for beam scanning will be 
deferred to a later report, but a typical circuit is presented as Figure 9. This 
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ANTENNA PORTS 



BEAN PORTS 


Figure 9-Siraultaneous beam-forming matrix for 16 x 1 (linear) array composed of 
coupler-type magic T elements and fixed phase shifters. 


particular configuration forms 16 independent orthogonal beams, individually or 
simultaneously, and the associated scattering matrix 16 jS £ . is given as Eq. (74). 
Ordinarily a beam-forming matrix of this type is associated with linear antenna 
arrays. 



A 

i 

i 

B 

1 

i j 

i 

‘ j 

4 

C. 

■ r ~ 

i 

D 


ij 

i 

1 j 


( 74 ) 


36 



S. . = e ^ in real time. Only the angle (+<//) is entered as an element of the 
scattering matrix here, for compactness of writing and is expressed in units of 
77 / 77 / 16 SlflCG 0*2 ~~~ U 77 ./PI 2 where ue (z odd). In Eq. (74), the partitions A. 
and D. . of matrix lg jS. . are null, representing matched conditions at all ports 
and perfect isolation between input or output ports. Since J6 jS. . is symmetric 
about the principal diagonal, B. . = C. . . These 16 x 16 matrices are partitions 
of the 32 x 32 matrix 16 jS. . . Matrix C is presented in complete detail in 
Figure 10. It is noted that 16 jS is a unitary matrix. Since( lfi S*.)t( S = I 

* J 10,1 ij y v 16,l ij) 



Figure 10-Partition C of matrix ,S. . 

n ij 16,1 ij 
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In like manner the scattering matrices of planar arrays can be written. If 


S. . is utilized as the basic component (matrix) , the composite beam-forming 

16,1 i j 

circuit of Figure 11 results; the associated scattering matrix being 


n 


2 * 



16,16 


s. . 


i j 


16 


A. . 


i j 


C. . 

i j 


k 
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Figure 11— Composite beam-forming circuit utilizing thirty-two ^ matrices. 



where the A. . = D. . =0 as before. Here B. and C. are not the same as in 
Eq. (74) since now each of these partitions is a 256 x 256 matrix of the overall 
512 x 512 matrix 16>16 S. ^ . Any column of C. . represents the phase conditions 
on the aperture of a 256-element square array of antenna elements. Only C. . 
is required to investigate various independent antenna beams, and B. . can 
always be generated by the simple expedient of transposing the subscripts (ij) 
in C. . . That is, B. . = C. . in general. If the antenna ports are used as beam 
ports the column matrices of B. j determine the phase gradients across the 
antenna aperture. Equation (76) illustrates the manner in which the aperture 
gradients can be identified for the planar array of 256 radiators depicted in 
Figure 12. The space relationship among the radiators is not indicated here. 
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(ANTENNA PORTS) 
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Number of internal phase shifters = 



Figure 12— Schematic representation of scattering circuit for 
256-element planar array of radiators. 


Under the transposition = C ., the column matrix B. 257 corresponds to 
the aperture phase gradient when the antenna port 257 is used as a beam port, 
for example, and 
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(77) 


The scattering matrix 1616 S. . is symmetric since it is equal to its own trans- 
pose, but the antenna ports cannot be employed as beam ports. Nevertheless the 
I6,i6^ij ma t r i x is reciprocal in the sense that either transmitting or receiving 
equipments can be utilized at the beam ports; a unique pattern exists for these 
two modes of operation. It is noted that the scattering matrices of simple coupler 
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and ring hybrids are both symmetric and unitary, as are the composite scatter- 
ing matrices ,S. . and ,,S... 

° 16,1 ij 16,16 ij 


BEAM COMBINING 


The possibility of weighted addition and subtraction of two or more of the 
independently formed beams derived from a planar aperture is now considered. 
Using the notation employed in the preceding analysis, the particular case of a 
16 x 16 element planar array is selected as an example and the composite beam- 
forming matrix of Figure 11 is utilized. Some attention must now be given to the 
absolute time delays involved in passing through the beam-forming matrix before 
performing any sum or difference operations. The beams are summed in space 
according to the superposition theorem discussed previously, therefore, Eq. (24) 
is repeated here for convenience, 


E = - l60I ° F ° — e _i ^ e"< VW f f f e ikR - ia;t . 
R 123 


(24) 


The surface on which the fields are observed will be taken as the sphere of 
radius R ; all source elements are assumed to be identical-, and the frequency of 
two simultaneously generated beams, A and B, is f = f = a: 2 rr . 


The interference phenomena on the sphere of observation now depend on the 
interaction of two beams 


and 


, - 1 ( $ j a+ ^2 A + S 3A ^ 

K A e 

. “d S i B + S 2B + S 3B ) 

K n t- 


f f f 

1 lA *2A *3A 

f f f 

1 IB l 2B l 3B » 


(78) 

(79) 


where k A and k B are constants of proportionality and equal or unequal weight 
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can be ascribed to these constants. By Eq. (23), the E A and E fi above can be 
rewritten as 


E =k r ~ i(S lA + S 2A ,| ’ S 3A ) sinn i^lA /2 sinn 2 y 2 A /2 Sint W 2 

A A sin ^lA /2 s i n 7 2 a /2 sin ^3A/ 2 


-i[(nj-l)y 1A / 2 + (n 2 -l )r2A / 2 + (n3 -l)y 3A / 2 ] (80) 


^B ~ k B 


-i($ 


1B + S 2 B + S 3B ) 


s io OjTj^/ 2 sinn 2 7 2B /2 
sin y 1B /2 siny 2B /2 


Sinn 3 y 3 B /2 

siny 3 B /2 


e ‘^ n l " 1 ) y \B /2 + ( n 2 ) r 2 B /2 + ) r 3B /2 ^ 


(81) 


For any point P(0 , 4>), and for beams formed by circuits of the type shown in 
Figure 11, the arguments of the exponential terms of Eq. (80) and Eq. (81) become 


1A 


+ S 2 A + ^3A 


(n 1 - 1) 

+ (ka, 


sin d cos d- 


+ & 


1A 


) +■ 


(n, - 1) 


■ (k a 2 sin 


6 s in <£ + a 2A ) 


(n 3 -l) 

+ 2 ( ka 3 cos 6 +a 3A ) (82) 
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(k Sj sin 6 cos g6 + a JB ) + 


(k a 2 sin 6 sin <p + a ?B ) 


(n x -1) 


IB 


2B 


3B 


(n 2 -l) 

2 


( n 3 - !) 

+ ^ (k a 3 cos Q + a 3B ) (83) 


A square array implies n 2 = n 3 = n ; a perfect image plane implies A “ ^ IB’ 
a iA = a iB ~ 7T ’ n i = Since all space terms cancel, the phase difference for the 
resultant far-fields of beams A and B, referenced to an end source element can 
be written as 


0 


= ? 


S 2B “ S 2A + S 3B 


S 3A + 


(n - 1) 


( a 2B “ a 2A + a 3B 


- a 3A ) 


(84) 


Referring to Figure 11, arbitrarily assume the selection of two adjacent 
beams so that a.. £ a„ n , but a, A a,_ . It is clear from the circuit that this 
implies that only one horizontal 16 t S ;j is energized, but all vertical 16 jS. 
are energized (transmitting case) and all radiators are energized. Also, since 
the adjacent beam signals enter the vertical 16 jS. at the same level, the delay 

S,_. The phase difference in this instance reduces to 

^ = S 2B - S 2A + ( a 2B ~ a 2A) • < 85 > 
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The term i p, above, is graphically presented in Figure 13 for 16 jS. when the 
adjacent beams are generated by columns 248 and 256 of C. . in 16 1( .S. . given 
in Eq. (76). The origin and physical significance of the terms (S 2B - S 2A ) and 
- n ~ — (a 2B - a 2A ) are based on the fact that aperture phase is referenced to the 
source element at the origin of coordinates. Although \p was derived from beam 
considerations in space, it is noted that all space factors are absent in the ex- 
pression for \p and only circuit parameters of 16 jS.. remain. It can therefore 



Figure 

be 


13— Graphic phase-center determination of 
ams using j S - - scattering matrix gradients. 
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be taken as the difference between the phase centers C A and Cg of the array 
beams A and B since 0 = C B - C A . Thus it can be seen that delay is a problem. 
If these beams are to be used to form a simple sum beam, a phase delay of (0) 
radians is required in the feeder of beam A . Similarly, (i p +v) radians of phase 
delay is required for a simple difference beam. 


APERTURE ILLUMINATION 

Originally uniform aperture-element illumination was assumed for each in- 
dependent transmitting beam formed by a simultaneous beam-forming matrix. 
This assumption was fundamental to the summation of the polygon of electric 
field components from the aperture, a geometric series. It is of some interest 
to investigate the aperture current distribution associated with the formation of 
sum beams as this provides intermediate beams in space and is prerequisite to 
more advanced concepts such as phased array monopulse. If the summation in- 
volves equal beam weights (i.e., equal principal maxima), equal coefficients can 
be assigned to the column matrices of C. . for each beam. Continuing with the 
planar example, the aperture current is obtained for the particular case of 
a 3 A “ a 3 B = 77/716 * a 2 A = 77 /16, a 2B r 3tt 16. For this case 

* = S 2B - 8 2 A + («2B - “ 2 A> < 85 > 

so that 


I 15 77 1777 (16 - 1) / 377 77 \ 

r 248 , 256 16 16 2 V 16 16/ 


The current amplitude distribution on the aperture is then obtained by summing 
the column matrices of C.. in [P kk ] [ 16 _ 16 S. J [P kk ] and [ 1616 S..] corre- 
sponding to the beams A and B, respectively, where P kk is a diagonal matrix 
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with all null elements excepting the principal diagonal. All P kk 1 except 
P = e j ^248,256 which corresponds to a translation of the reference 

256 ,256 

plane in forming the A beam by excitation of input (256) of 16 l6 S . . . The B 
beam is formed by direct excitation of input (248). 

Each element of the column matrix is therefore delayed in phase by an 
amount ^ 248 256 = 137 t/ 16 . The current amplitude (I) on the aperture when 
beams A and B coexist, and are properly phased, is given by a sum of column 
matrices: 


C 257,248 
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' ^2 48,25 6 
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*258 
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5 11,256 

^248,256 


*511 
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512,256 

“j ^24 8,25 6 


*512 

— — 


— 

— 
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Figure 14 illustrates the discrete transmitting aperture current amplitude dis- 
tribution for the sum of beams A and B generated by the previously assumed 
values of a, and a . The current is uniform in the a direction, but has a 
cosinusoidal envelope in the a 2 direction. This result holds for adjacent beams 
with equal excitation. The current amplitude can be plotted directly, now, by 
generating cross-sections of the distribution using the simple rule that the angle 
between source element current of beams A and B is given by 

- x (a B - a A ) (87) 

where a represents the slope of the linear aperture phase plot of Figure 13, and 
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x is the distance to a radiator from the array center in units of a 2 . Cuts of 
Figure 14 also represent the current distribution for a linear array of 16 
source elements for equal-weight beam addition. The corresponding sum and 
difference electric field patterns associated with this current distribution are 
shown in Figure 15. 

Certain interesting similarities and differences exist between the transmit- 
ting and receiving modes of operation. In forming a single transmitting beam 
with the 16 x 16 planar array in the previous example, the aperture illumination 
is uniform and a single input energizes 256 output terminals of s ... In 



2 



a 2A’ a 3A'a- 3B - -j| RADIANS 
a 2B * RADIANS 


Figure 14— Envelope of planar aperture current 
distribution resulting from summation of a pair 
of equally weighted, adjacent, orthogonal 
transmitting beams. 
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Figure 15— Sum and difference pattern formation of 16 X 1 (linear) array of half-wave 
dipole elements (H-plane cut), a j - a 2 ~ A./2. 

receiving a single planar wavefront there exist 256 equal signals at the aperture 
and (256) signals at the beam ports. The transmitting and receiving patterns are 
identical for these two cases, as are the aperture current distributions. It is 
noted that the current distribution on the aperture is not necessarily the same 
for all transmitting/receiving modes for which the patterns are identical. A 
first-order cosine beam is associated with uniform aperture illumination on 
reception, but cosine-uniform illumination (as in Figure 14) during transmission 
if two adjacent beams are summed with equal weight in a planar array. 
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CONCLUSION 


This report represents an attempt to obtain a realistic description of the 
far-field radiation patterns for a uniform planar array and, in particular, con- 
cerns itself with combinations of gradients obtained from a passive multiple 
beam-forming matrix. Although the array factors mentioned herein can be 
found in many textbooks, and elsewhere in the literature, the viewpoint that a 
uniform planar array is representable as a three-dimensional lattice of oscil- 
lators appears to be less common in the antenna literature. For many purposes 
this is a useful notion. It is possible to analyze the position of beam maxima, 
nulls, crossover levels, etc. without actually plotting the field patterns. Distinc- 
tions are easily made between non- realizable beams (or spurious beams) and 
grating lobes. The role of the element source pattern can be seen very clearly 
and the artifice of isotropic element source patterns can be avoided. Finally, 
the diffraction power-patterns are formed for physically realizable element 
source factors and three lattice factors, and are then integrated over the hypo- 
thetical unit sphere to obtain the total radiated power for various combinations 
of element spacing and phase gradients. A weighting or normalization of the 
diffraction patterns is then applied consistent with energy conservation, pro- 
viding a realistic measure of directive gain. 

The introduction of the scattering matrix notation for multiple-beam for- 
mation by passive circuits is useful to display the natural order in the signal- 
processing circuit. Scattering concepts are of proven value in many physical 
disciplines and particularly in the analysis of electrical circuits. They are not 
essential to array analysis or to the understanding of multiple beam-forming 
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matrices. Although the full benefit of the scattering notation may not be evident 
from the limited discussion given in this report, an examination of the cited 
references makes clear the extent to which the basic concepts can be carried 
in the analysis problem. As a specific example scattering matrices can readily 
be converted into transfer or chain matrices which admit the possibility of 
tandem-connected circuits. Many other advantages could be cited. These lead 
to a specialized discussion of eigenvalues, eigenvectors, canonical forms, Pauli 
spin operators, and related topics in linear algebra and topology. The value of 
formulating the beam-forming technique is that a very large body of mathematical 
operations becomes applicable. 

This report is planned as the first of a series of reports on the subject of 
phased arrays and passive scattering matrices, and represents an attempt at 
unifying some of the topics of elementary planar antenna array system analysis. 

It is intended to serve as a point of departure into more sophisticated treatments 
of the overall subject. Less idealized arrays can be investigated within the 
framework of the notation developed, and more complex systems can subse- 
quently be described by utilizing the basic concepts presented. Subjects such 
as phased array monopulse, high-order cosine taper illumination, aperture 
gain, reactive aperture power, array energy conservation, and mutual coupling 
phenomena constitute logical extensions of the present report. An extensive 
study of microwave scattering junction topology and basic coupling techniques 
appears essential to achieving increased design flexibility in phased array 
systems. 
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Appendix A 

Radiation Pattern of Crossed-Dipole Antennas 


The radiated fields of the crossed dipole antenna shown in Figure (1) can be 
determined by the Poynting integration method. 1 The radiation vector from a 
linear antenna is given by 


where 17' 


N = 



* dl 


( 1 ) 


radius from origin to differential element, d'C, on the antenna 


7 1 = radius from origin to the point P at which the field is to be 
calculated 


0 = angle between r and 7' 


and the time dependent factor e 1 " 4 is understood (i.e., omitted). The vector 
potential at point P is 


M ■ 


kr 

477 r 


N 


( 2 ) 


and the fields at point P are given by 


E = V(V-A) - j a A 
k 2 


B = V x A 


(3) 


^Ramo, S. and Whinnery, J., "Fields and Waves in Modern Radio”, New York: John Wiley and Sons, 

Inc., 1953- P- 505 ff. 
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The only components of the above fields which do not decrease faster than 1/r 
are 


E* = 


Zl2± i c-ikr 
477 r 


N 0 



477 r 


e”'> kr 


H, 


= JlL 

477r 


‘ jkr N, 


H. 


ZjJi e -jkr 

47Tr 


N, 


For a specified antenna configuration and current distribution, equations 
(1) and (4) can be solved for the components of the radiated fields. For the 
crossed-dipole antenna with half-wave dipole elements and an assumed sinus- 
oidal current distribution, the radiation functions (components of the radiation 
vector) are shown in Appendix A to be 


2I o 

N «(# a ,#,#> = — 

a k 



(sin 0 cos 0 sin 


-sin#cos# a ) 


+ 




(s i n <fi cos # co s 8 g + sin#sin# a ) 


2I o 

tJ > <£) - — COS 4> 


sin# cos 


(f) 


CO s d cos 




1 - a 2 


1 -/3 2 


where y = difference in phase of the currents on the individual dipoles 


a = cos#cos# + sin# sin# sin 0 


/3 = cos # s in # fl - s in # cos # a sin 0 


# = tilt angle defined in Figure 1 





Substituting for the tilt angle and phase difference in Equations (5) and (6) 
yields the components of the radiation vector for the desired antenna orientations 
and polarizations: 


1) 6 =0°; linear polarization (y = 0°) 


2I o 


77 

COS ( — cos 


S 1 


in 9 


sin <£> cos 9 cos ^-sini9sin<£ 
1 - sin 2 9 sin 2 0 


(7) 


2 j cos 0 cos - sin 0 sin 

N<*(0> 6 , 4 >) - r— V2 


77 


K 1 - s in 2 t? s in 2 0 

The substitution of 4 > = 0° in Equations (7) and (8) leads to the radiation 
functions in the E -plane (0=0° plane) as 


(8) 


N 5 (0, 0 , 0) r~ 


r, T cos I — cos 6 
2I o \ 2 


s 1 


n & 


(9) 


2I 0 

N 0 (O, 0 , 0 ) = -— 2 


( 10 ) 


The substitution of 6 = 90° in (7) and (8) leads to the radiation functions in 
the H-plane (8 = 90° plane) as 

2I„ 


Nn (0, 90 , 4 >) 


( 11 ) 


(0, 90 , 4 >) = - — 


2I o COS (i Sin 


COS 4 > 


(12) 
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$ a - 0°; circular polarization (y = 90°) 


21 , 


N.(O,0,0) = - — - 


cos[— cos s in </> cos 6 cost — s i n 6 s in <b 

' 2 / . \ 2 j 

+ J — 


\ 


in 6 


sin 


1-sin 2 6 sin 2 


21. 


N ^(°. 9, 4>) = j — — cos 4> 


77 


cos ( — sin 6 sin cb 
1 2 


1 - sin 2 6 sin 2 4> 


The radiation functions in the E-plane are 


N*(O,0, 0) = - — 


OT cost — COS 6 
2i o \2 


77 


s i n 6 


VO, e, o) = j 


21, 


The radiation functions in the H-plane are 


N e (0, 90°, 0) = - 


21 , 


N^(O,9O°,0) = j 


or cos — sin d: 
Zl o \ 2 


3) # a = 45°; linear polarization (y = 0) 
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( 20 ) 


The radiation functions in the E -plane are 


( 45 °, #, 0 ) =, - 


2 v2 L 


sin # cos 


(it cos ") 


1 - — cos 2 6 
2 


(21) 


N, (45°, 6, 0) = 0 


( 22 ) 


and the radiation pattern is plotted in Figure 4. 
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The radiation functions in the H-plane are 


N g (45°, 90°, 4 ) =- 


2 V2 


COS 


7 T 


\2V2 


sin cp\ 


0 1 

1 sin 2 4 
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(23) 


N ^(45°, 90°, 4 ) = 0 


(24) 


and the radiation pattern is plotted in Figure 5. 


4) 6 = 45°; circular polarization (y —n /2) 
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The radiation functions in the E -plane are 


N,(45°,£,0) = -lli£ (l.j) 


sin 0 cos 


{jTT cos e ) 


1 - i cos 2 9 
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(27) 


N # (45”,«.0)=^io (I + j 
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cos { =- cos 6 
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(28) 


and the radiation patterns are plotted in Figure 6. 


The radiation functions in the H-plane are 


/2 1 o 
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and the radiation patterns are plotted in Figure 7. 
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Addendum to Appendix A 

Derivation of Radiation Vector for Crossed-Dipole Antennas 


For the k/2 dipole oriented along the Z'-axis (Figure 2) the radiation 
vector is given by 


N,=Jl 


i \ e i k l 7 '|cos0 


J k/4 

V e 

-k/4 


j k I z * I cos </» dz » 


If u and a 2 , are unit vectors along r and 7 ' , respectively: 


u = a x sin 0 cos 0 + a y sin#sin0 + a z cos ^ 
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z y a z a 
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For a sinusoidal distribution of current 
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If the current for the X./2 dipole along the Z"-axis (Figure 3) is assumed to be 


i z „ - a z „ I Q e“ ir sin 


then its radiation vector is determined in like manner to be 


2I n e"j y /„ \ 

cost — /3) [-a cos 0 + a sin 0 

k(l-/J») \ 2 ) y 


The total radiation vector is 
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and the spherical components are: 
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Figure 6-Radiation Patterns on E-Plane for t J A =45°, y - 90 ° . 



Figure 7-Radiation Patterns on H-Plane for 0 A = 45°, > = 90°. 
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